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In this paper we study the physics potential of the KamLAND detector in probing neutrino 
oscillation parameters through observation of supernova neutrinos. In particular, we discuss the 
fT^ ' possibilities of probing the mixing angle 613 and determining the sign of Am32 from the total 

' charged current(CC) event rates on the proton and ^^C target, as well as from the CC spectra. We 

. discuss the chances of probing the earth matter effect induced modulations from the observation 

04 ■ of CC spectra in the different CC reactions in KamLAND and find the volume required to get a 

^ ' statistically significant signature of the earth matter effect in different energy bins. We also calculate 

^ , the event rates expected in the neutral current (NC) reactions on Carbon and free proton and 

■ investigate if the charged current to neutral current ratios, which are free of the absolute luminosity 

' uncertainty in the supernova neutrino fluxes, can be useful in probing the oscillation parameters. 
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CN . PACS numbers: 14.60.Pq, 13.15.+g 
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I. INTRODUCTION 



' KamLAND is a 1 kton liquid scintillator detector situated in the Kamioka mine in Japan. This detector is at 
OA ', present engaged in measuring the flux of electron antineutrinos coming from the various nuclear reactors in Japan. 
^""^ ' In 2002 it reported a phenomenal result H^, providing the clinching evidence for the Large-Mixing- Angle(LMA) 
\ MSW 0] solution to the solar neutrino problem - which was emerging as the favoured solution to explain the 
. solar neutrino shortfall in Homestake, Kamiokande, SAGE, GALLEX/GNO 0|, Super-Kamiokande(SK) 0j and the 
' Sudbury Neutrino Observatory(SNO) H 0, experiments. The combined two-neutrino oscillation analyses of the 
Qh, solar neutrino and KamLAND data including the most recent SNO salt results give the best-fit values of neutrino mass 
O ^ squared difference and mixing angle as Atoq(= Am|i) = 7.2 x 10^^ eV^, sin^ Oq{= s\v? 812) = 0.3 j2ll3|. The neutral 
' current data from SNO had been instrumental in ruling out the maximal mixing and dark-side solutions (0f7) > 7r/4), 
I implying Arri^i > a result now confirmed at more than 5a level by the latest results from SNO [3,|a|- 

^ I Compelling evidence in favour of neutrino oscillation have also come from the observation of up-down asymmetry 
• , in the Zenith angle distribution of atmospheric neutrinos in the SK detector. The best description to the data comes 
' in terms of Vr {v^ Vt) oscillations with maximal mixing and Am^(j„(= ATO32) = 2.0 x lO"'^ eV^ from a 

recent re-analysis of data performed by the SK collaboration [lH . 

In the realistic three neutrino oscillation scenario the two sectors - solar and atmospheric - get coupled by the 
mixing angle ^13, which is at present constrained by the reactor data as sin^ 6'i3 < 0.03 at 90% C.L. This is 

only an upper bound and one of the major issues to be settled by future experiments is - how small is ^13. Another 
important question which has still not been ascertained by the current data is the sign of Am§]^. Both these issues 
are expected to be addressed in future long baseline experiments 0, Q| . Reactor experiments have also been shown 
mising capat 
core collaps 

the above two issues 0, |3 ^3 ■ The feasibility of supernova neutrino detection has already been demonstrated by 
the chance recording of the neutrino events from the supernova SN1987A by Kamiokande [2l| and 1MB [23| detectors. 
The few neutrino events detected were sufficient to confirm the basic theory of core collapse supernova and at the 
same time could be used to set limits on neutrino properties. SN1987A occurred in the Large Magellanic cloud at 
a distance of about 50 kilo parsec(kpc) from earth. For a supernova in our galactic center, we can assume a typical 
distance of 10 kpc and therefore expect many more neutrino events. In type II supernovae the explosion takes a few 
hours to penetrate through the envelope and only after that the optical and other electromagnetic signals emerge. 



experime 

recently to ha ve p romising capabilities to probe ^13 [Tsl IT^ and under very specific demanding conditions also the 

^32 



sign of Ato^j jl7j | . A core collapse supernova (SN) in our galaxy can also play a significant role in throwing light on 



Whereas neutrinos start coming out right after core bounce and are therefore useful for alerting supernova telescopes. 
Therefore it is important that the present detectors be ready and the potentials fully explored before the next galactic 
supernova occurs. The capabilities of present and future detectors towards detection of neutrinos from a galactic 
supernova has been investigated by many authors (see e.g. for a list of references). In particular, the ability of 
SN neutrinos to probe oscillation parameters have been considered in many papers [isL IT^ l20l l2 ^ . For example, it 
is shown in (18.. 19, 20J that if the propagation of neutrinos is non-adiabatic in the SN matter, then they provide 
us a window to glean into the important parameter ^13. Since both neutrino and antineutrino fluxes arrive from 
a SN, one can in principle also get an idea about the neutrino mass hierarchy, by observing the difference in the 
neutrino and antineutrino event rates in the detectors. The possibility that the neutrinos can in addition undergo 
earth matter effect depending on the location of the SN opens up interes ting possibilities of studying earth matter 
effects by comparing the neutrino fluxes at different detectors [2^ l28| . This in turn can help to distinguish 

between the mass hierarchies |0, . In a realistic picture one should consider the uncertainties coming from 

the SN neutrino fluxes as discussed in detail in . However one can construct variables through which one can still 
hope to discover 2 — 3ct effects ^ . 

In this paper we do an extensive study of the KamLAND detector towards detection of SN neutrinos. The dominant 
reaction for detection of neutrinos in KamLAND is the i>e + p + n. The liquid scintillation detector material 

also offers the possibility of studying the neutrino interactions on the Carbon nuclei |^ . This can proceed via both 
charged current - ^^C(z^e, e)^^iV and ^^C(i/e, e+)^^i? - and neutral current - ^^C{iyx,'^xy^C* . The charged current 
reactions on ^^C are detected by the delayed coincidence of the final state positron/electron with the electron/positron 
produced by subsequent /3 decay of the daughter nuclei. The event signatures are separated by about tens of milliscc 
in time and can be tagged. The final state carbon produced in the neutral current reaction on ^^C produces a 
15.1 MeV 7 - the detection of which requires large homogeneous detector volume and good energy resolution [s^ . 
KamLAND detector satisfies both these criteria. The low energy threshold of KamLAND also makes it possible to 
detect the neutrino-proton elastic scattering reaction for the SN neutrinos (s^. Neutrinos and antineutrinos of all 
three flavors participate in this neutral current reaction and the event rate is comparable and can even be greater 
than the i>e + p events. 

We explore the information that one can obtain on neutrino oscillation parameters from the observation of supernova 
neutrinos through the above mentioned reactions in KamLAND. Since the average energies of the v^^lvTiy^ilvT) are 
larger than the average energies of i/e(Pe), flavor oscillations cause enrichment of the ^'e(S'e) flux in higher energy 
neutrinos, thereby resulting in the hardening of the corresponding spectra. Therefore, one expects an enhancement in 
the CC events in presence of oscillation. This enhancement could be particularly prominent for the CC reactions on 
^^C, for which the energy thresholds are very high and for the unoscillated v^. or v^. one expects very few events. Since 
the degree of enhancement depends on the survival probabilities and hence on the oscillation parameters, especially 
the mixing angle ^13 and the neutrino mass hierarchy, the total observed event rates in KamLAND can be used to 
constrain these parameters. We study the variation of the total event rates with ^13 and discuss the limits on ^13 
that one can obtain for a given mass hierarchy. We also expound the possibility of differentiating between the two 
hierarchies from the observation of energy integrated event rates, thereby determining the sign of Am|]^. 

The hardening of the z^e(f'e) flux is also manifested in the observed electron/positron spectrum in the detector. 
Thus, by analysing the observed CC spectrum, one may gain information on oscillation parameters. In addition, if 
the SN neutrino cross the earth matter to reach the detector, then the survival probabilities would be modified due 
to matter effects inside the earth. This "earth regeneration" effect is dependent both on the solar neutrino oscillation 
parameters, as well as on 6*13 and the hierarchy. Thus the observation of these earth modulations would also give 
insight into these parameters. The earth matter effect is seen as very fast oscillations superimposed on the slowly 
varying spectrum of the neutrinos coming from the SN. Since the variation with energy is fast, averaging effect due 
to energy resolution of the detector could wash out the earth effect completely from the energy spectrum. Therefore 
one needs a detector with very good energy resolution. KamLAND being a scintillation detector, has a very good 
energy resolution. We compare the degree of earth matter effect in the various CC reaction channels in KamLAND . 
The reactions involving the ^^C are seen to have large earth effect induced fluctuations in the lower energy bins. 
However, with the present volume of 1 kton, KamLAND lacks the statistics for drawing precise constraints from the 
^^C spectral data. We present the minimum detector mass required to detect these modulations at a statistically 
signiflcant level in various energy bins. We also present the spectrum for a 50 Kton scintillator detector and discuss 
the chances of probing ^13 and hierarchy. 

Finally, we discuss the role of the uncertainties originating from our lack of knowledge of supernova neutrino fluxes. 



An alternative proposal was discussed in 1291 1301 , whereby using the known values of the neutrino oscillation parameters one could 
determine the SN parameters. 



The NC reactions (both on Carbon and free proton) give a measure of the total neutrino flux coming from SN. We 
introduce the the ratios of charged and neutral current events in which the absolute normalisation uncertainties, e.g. 
due to the SN absolute luminosity or the distance of the supernova from earth, get canceled. We discuss the merits 
of these ratios in studying the oscillation parameters. We find that these ratios provide a good handle to probe 
oscillation parameters particularly when we consider the neutral current v — p scattering reactions. 



A star in an advanced stage of burning has an onion shell structure with different nuclear reactions going on in 
each shell. After millions of years of burning the nuclear reactions stop inside the innermost shell (i.e. the core) with 
matter consisting mostly of ^^Fe like nuclei. Once this energy generation in the core stops and if the core has mass 
greater than the Chandrasekhar mass, the star undergoes a very fast gravitational collapse. During this stage, electron 
neutrinos are produced through neutronisation reactions - electron capture on free protons and nuclei. During the 
early stage of the collapse, matter density of the core remains below the "neutrino trapping density" ( ^ 10^^ g/cc), 
the neutrinos have mean free path much larger than the radius of the core and hence escape. Within a duration of 
tens of milliseconds, this collapse compresses the inner core of the star beyond nuclear matter densities. Consequently 
this inner core rebounds and creates a shock wave deep inside. This shock wave propagates radially outward through 
the star, gets stalled in the middle due to loss of energy through nuclear dissociation, eventually gets re-energized 
and finally hits the outer mantle in a few seconds, supplying an explosion energy of a few times 10*^^ ergs. This is 
believed to be the cause of the type II supernova explosion. During the shock phase, due to high temperature inside 
the core of the star, thermal v and v are produced in all three flavors from e~^e~ pair production and "urea processes" 
[s^ l . This production of neutrino pairs stops when neutrino density becomes high enough so that the inverse process 
of pair production balances the direct process. The total energy of the neutrinos from pair production at the shock 
phase is much larger than that of the neutronisation neutrinos. Almost all the gravitational binding energy released 
due to the collapse is radiated away through these neutrinos with only a few percent going to the explosion. These 
neutrinos inside the core of the star are in equilibrium with the ambient matter density and their energy distributions 
are close to Fermi-Dirac, with the characteristic temperature given by the temperature of last scattering. This is 
seen through simulations and through the analysis of SN1987A neutrinos. Out of all neutrinos, f^, Vt and their 
antiparticles interact with matter only through neutral current, whereas z/g and 9^ have both charged current and 
neutral current interactions. Since matter inside the SN is neutron rich, i/g's interact more than DeS. Therefore v^,, 
Vr and their antiparticles decouple first and hence have the largest temperature, followed by the Ve, while Ve decouple 
last and have the smallest temperature. In this paper we consider the detection of neutrinos arriving at earth from the 
cooling phase of the SN burst, and assume that all the three types of neutrino gas {ve, Vf. and Vx, where x stands for 
the four species v^, i/^, Vt and i>r) have the Fermi-Dirac energy distributions with equilibrium temperatures as 11, 16 
and 25 MeV respectively (see for instance js^J). However although the temperature hierarchy are generally believed 
to be in the above order the exact values of the average energies may differ depending on the details of supernova 
simulations jlE 1^ • Such simulations also indicate that the high energy tail of the neutrino spectra from SN is 
better reproduced by a "pinched" Fermi-Dirac distribution with a non-zero chemical potential [39j. Finally in this 
paper we have considered that the total binding energy (3x 10^^ ergs) released in the core collapse of the star is equally 
partitioned between various flavors. But detail simulations may give rise to scenarios where the v^, and Uf, luminosities 
can vary between (0.5-2) times the Vx luminosity 37]. The propagation of the neutrinos and their conversion from 
one flavor to another depend on matter density proflle in the supernova mantle. A radial profile of 



Neutrinos and antineutrinos of all flavors are produced in a supernova burst. Before reaching the detectors, the 
neutrinos travel through the supernova matter, propagate through vacuum and finally may even travel through the 
earth's matter depending on the time of the SN burst. For the range of neutrino oscillation parameters consistent 
with the solar, KamLAND reactor and SK atmospheric data, one expects large matter enhanced oscillations inside 
the SN and regeneration effects inside the earth. The neutrino detectors on earth can detect individual fluxes of Vf. 
and Vf, through the CC reactions and the combined fluxes of all the six neutrino species through the neutral current 



II. NEUTRINO EMISSION FROM TYPE-II SUPERNOVA 




(1) 



with C~ 1-15 provides a good description of the matter distribution for p > 1 g/cc psl liollZll li ^ . 



III. CALCULATION OF THE SURVIVAL PROBABILITY 



reactions. For the case where one has only active neutrinos, it can be shown that the probabihties involved in the 
calculation of the event rates in the earth bound detectors are only the survival probabilities Pee and Pgg. 

In the three-generation scheme ^ neutrinos would encounter two resonances in the supernova matter. The first one 
corresponds to the higher atmospheric mass scale at a higher density and the second one corresponds to the lower solar 
mass scale at a lower matter density. The higher resonance is controlled by Am§2 and ^13 while the lower one is mainly 
controlled by /S.m\i and 612. The transition probability depends crucially on the sign of the mass squared difference. 
Therefore depending on the sign of Am32, which is still completely unknown, we can have two possibilities for the 
neutrino mass spectrum: Normal Hierarchy (NH) when ATO32 > 0, or Inverted Hierarchy (IH) when Am|2 < 0. We 
briefly summarize below the expressions for the survival probabilities in the three-generation scheme of the neutrino 
mass spectrum, with the current allowed mass and mixing parameters and for both normal and inverted hierarchy. 
We use the standard Pontecorvo-Maki-Nakagawa-Sakata parametrisation for the mixing matrix [45j |. 

A. Normal Mass Hierarchy (NH) 

The mixing angles for neutrinos in matter for normal hierarchy are given as 

tan2 2rVrl - Am^^ sin 2^12 



i2 26l7J('"^ - Ato2^ sin 2013 
A{r) being the matter induced term given by 



tan2 20r3M = 2 ^ (3) 

' Am§i cos 26*13 - A(r) 



A{r) ^ 2V2GFNAne{r)E, (4) 

where Gp is the Fermi constant, Na the Avogadro number, ne{r) the ambient electron density in the supernova at 
radius r and E,y the energy of the neutrino beam. Since A >> Am3i >> Ato|i at the point of production of the 
neutrinos, from equations ^ and Q we have, 9^ « -I « O"^, which means that neutrinos are created in pure i'™ 
states and the expression for the survival probability is 

Pee = PnPLP^e + Ph{1 - Pl)PI + (1 - Ph)PZ (5) 

where P® are the Vi — )■ Ve transition probabilities and may depend on earth matter effect if neutrino cross the earth. 
Ph and Pl are the "jump" probabilities at the high and low density resonance transitions respectively and for these 
we use the standard double exponential forms ^Jj, ^4^] . Although this form derived in was in the context of solar 
neutrinos and exponential density profile, it was shown in that this holds good for generic SN density profiles as 
well. 

For antineutrinos the expressions (01 and Q hold with the minus sign in the denominators replaced by a plus 
sign. Then, at the antineutrino production point 0^ « w (i.e. they are created in the state ) and the anti 
neutrinos the survival probabilities are given by 

P-e-e = {l-PL)P?e+PLPfe (6) 

where the jump probability Pj, for the anti neutrinos is given by the standard expression as given in [43. |49|. 



B. Inverted Mass Hierarchy (IH) 

For inverted mass hierarchy (Am|i ATO32 < 0) the 12 matter mixing angle is same as that for NH both for 
neutrinos and antineutrinos. The 13 mixing angle in matter is given by the same expression as for NH with the sign 



^ We disregard the positive evidence from LSND experiment Ksl . Study of SN neutrinos in the context of four neutrino species have been 
done in 1441 . We also neglect the CP violating phases which is a justified approximation as for equal luminosities and temperature of 
the /i and r neutrinos coming from SN, the relevant probability at the detector is only the electron neutrino survival probability which 
does not depend on the CP phases. 
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FIG. 1: Ph, Pl as a function of Am|i and Ph as a function of Am|i at neutrino energy — 20 MeV 



of Am|]^ reversed. Then at the production point of the neutrino ~ J and ^"3 ~ and are produced as and 
the neutrino survival probabihty is 



Pee = PlP^, + {1-Pl)P^ 



(7) 



whereas at the antineutrino production point ^"2 ~ and ^"2 w -I so that v,. are produced as i'™ states and the 
survival probability is given by 



Peg = {l-PL)PHPt + PLPHP^e + {^-PH)Pi 



3e 



(8) 



C. Dependence of survival probability on oscillation parameters and energy 



The dependence of the jump probabilities Pl and Pl on the 1-2 mass and mixing parameters and of Ph on the 1-3 
mass and mixing parameters has been shown in the figure ^ for an illustrative neutrino energy of 20 MeV. We see 
from this figure that for Am^i > 10~^ eV^, Pl = = Pl for all values of 6*12. We have checked that this conclusion 
remains valid for all energies of the SN neutrinos. Since KamLAND and the solar neutrino data have independently 
confirmed that Arn^i lies in the LMA region, we can safely take Pl = in the expressions for the probabilities. For 
the current range of Arri^i allowed by the atmospheric neutrino data, Ph sharply depends on the value of ^13. It is 
therefore useful to consider the survival probabilities in three limiting range of 6113 ^ 



^ In what follows, we keep sin'^ 612 fixed at its current best-fit value of 0.3, and Ar?igj fixed at 0.002, the best-fit value presented by the 
latest re-analysis of the SK atmospheric collaboration. The probabilities depend on Arrij^ if there are earth matter effects, and we keep 
Arriji fixed at its best-fit value of 7.2 X lO^^eV^. 




FIG. 2: Plot of Pee.^ , Pii^ , Pl^ and P//^ vs tan^ 613 for three different values of neutrino energy. The left-hand panels show 
the case for neutrinos coming straight from the SN while the right-hand panels correspond to the case when they travel earth's 
mantle and core as well. The other oscillation parameters are kept fixed at their current best-fit values. 



• The small 9iz range: For tan^ 6*13 ^ 10 ^, Ph = 1 practically over the whole relevant energy range of the SN 
neutrinos, except for very small energies. This is the regime of extreme non-adiabatic transitions. 

• The intermediate 613 range: In the range 10^^ ^ tan^ ^13 ^ 10^"^, Ph ranges between and 1, depending on 
the exact value of ^13 and energy. 

• The large 613 range: For tan^ 613 ^ 10~^ Ph — for all energies, in the entire range of Amgj^ allowed by the 
atmospheric neutrino data. This is the region of complete adiabatic propagation. 



1. No earth matter effect 

For this case 

P® = COs2 013 COS^ 012 (9) 
P® ^ cos2 013 sin2 012 (10) 
P3® = sin2 0i3 (11) 

Depending on the value of 0i3, the probabilities assume the following forms: 

• In the range tan^ 0i3 ^ 10^^, Ph = 1, sin^ 0i3 sa so that, 

PiI^ = sin2 0i2=P/f (12) 
Pi^^ = cos2 0i2=P/f (13) 



• In the range 10 ^ tan ^13 ^10 , Ph is non-zero but sin ^13 can be taken as zero and, 

P^/ = Ph sin2 ; Pi" - sin2 0i2 (14) 
Pg" = cos^ 012; P/f = Ph cos^ 0i2 (15) 



• For tan^ 0i3 ;^ 10"^ Ph = 0, sin^ 013 can be non-zero ■* and the probabilities are, 



Pg" = sin^ 013 ; P/f ^ cos^ 0i3 sin^ 0i2 (16) 
Pi^^ = cos^ 013 cos2 012 ; P/f = sin2 0i3 (17) 



The dependence of the probabihties on 0i3 for the case of no earth matter effects are shown in the left-hand panels 
of figure 12 for three different values of the neutrino energy. From the figure |21 and the equations for the probabilities 
at various limits we note the following points: 

1. In the small 0i3 regime (tan^ 0i3 ^ 10^^) the survival probabilities for both and are independent of the 
value of 013 and neutrino energy. The flavor conversion is governed only by the solar mixing angle 0i2, and one 
cannot distinguish between normal and inverted hierarchies in either the neutrino or the antineutrino channel. 

2. For IH(NH) the survival probabilities for Ve{i^e) do not depend on Ph- Hence they remain largely independent 
of the value of 0i3 and the neutrino energy and depend mainly on the solar mixing angle 0i2 in the entire 
oscillation parameter range considered in this paper. Dependence on 0i3 comes only for NH(IH) in the Vei^e) 
channel. 

3. In the intermediate 0i3 range (10~^ ^ tan^ 0i3 ^ 10"'^) and with NH(IH) as the true mass hierarchy, the 
neutrino (anti-neutrino) survival probabilities depend sharply on the value of 0i3 and neutrino energy, through 
the dependence of Ph on these variables. 

4. In the large 0i3 range (tan^ 0i3 ^ 10~^ ) and with NII(III) as the true mass hierarchy, the neutrino(anti-neutrino) 
survival probabilities are very small, within the CHOOZ limit of sin^ 0i3, and do not depend on energy. 

5. In the case of Ve{ve) for NII(III), the survival probabilities decrease from sin^ 0i2(cos^ 0i2) to ~ sin^ 0i3(sin^ 0i3), 
as tan^ 0i3 changes from ~ 10~^ to ~ 10~^. Since the best-fit value of the solar mixing angle is sin^ 0i2 ~ 0.3, 
Pgl^has a much sharper decrease with tan^ 0i3 than Pg^ ■ Therefore the Ve events could be expected to have 
better sensitivity to 0i3 on this count. 

6. In the range 10^^ ^ tan^ 0i3 ^ 10"'^, Pg^ and P/g^ depend on Am§]^ through Ph but over the currently 
allowed atmospheric range of this parameter, the dependence is not very strong. 

7. The survival probabilities in general depend on the value of 0i2, excepting Pg^ and Pg|^in the large 0i3 regime, 
but there is no dependence of the probabilities on Arn2i- 



2. Including the earth matter effect 



For this case 

P® = cos2 0i3K(25en)|2 (18) 

P® = cos2 0i3|A®(25en)p (19) 

P3® = sin2(?i3 (20) 

j4® (2gen) and A®g(2gen) are the Vi Vg, transition amplitudes inside earth for two generations 113,1^, with matter 
term A{r) replaced by yl(r) cos^ 0i3. Since the supernova neutrinos have typical energy of ^ 10 MeV the matter 
potential in the earth is << the high mass scale Am§]^. However, oscillations due to the Amji scale are affected by 
the earth matter effect. 



* However if we restrict ourselves up to the CHOOZ bound then sin^ 613 stays very small. 



• For tan^ 613 ^10 ^, Ph = 1, sin^ 9i3 can be set to zero and the probabilities are, 



= ^2®(2.9en) = P/f (21) 
FiI« = Pfe(25en)^P/f (22) 



• For 10 ^ ^ tan^ 6*13 10 3 the probabiUties can be written as, 

pNH ^ p pe . piH ^ pe (^o\ 

ee ~ ^H^2e i ^ee ~ 2e 

^'i^^-A® ; ^'/f^^^A® (24) 

where we have used sin^ ^13 ~ which is a good approximation in this range. 

• For tan2 6'i3 lO'^, sin^ 013 may be small but non-zero in the range allowed by the CHOOZ results. However 
as we have seen before that Ph — and the probabilities arc, 

Pil^ = sin2 0i3 ; Pif = P2® (25) 
P,V = P?e ; Pif-^Ois (26) 



We can define the following quantities valid for all values of 6113 : 

Rg" = Ph cos^ 013/,,,, ; Pff = cos2 ^13 f^eg (27) 



Rg" = -C0S2 013 /.eg ; i?|f = -PffCOS^ 013 freg (28) 



where Rf^ = Pj^(with earth effect) - P^^ (without earth effect), i can be e or e and x can be NH or IH and 

freg = - siu^ 0i2 (29) 

freg = " siu^ 0i2 (30) 

The freg and freg Carry the information about the regeneration effects inside the earth and are dependent on the 
neutrino energy, 0i2 and Am2i. 

From Eqs. p7(l and 128|l we see that while R^^ and Rg^ are independent of Ph, Ree^ and P||^ are linearly 
proportional to it. Thus, while the Ve(ye) spectrum will have an energy (and 0i2 and lS.m\^) dependent, but largely 
013 independent modulation induced by fregifreg) for NH(IH), for IH(NH) the degree of earth matter effects will 
depend sharply on the value of 0i3. In particular, when Ph — 1 the earth effects will be same for both the hierarchies. 
However, as 0i3 increases, Ph decreases, resulting in smaller earth regeneration effect in the S'e(i^e) spectrum for 
IH(NH). The earth effect completely vanishes in the Peii^e) channel for IH(NH) when Ph = for tan^ 0i3 ^ 10"'^. 

We show the impact of earth matter effect on the survival probabilities in the right-hand panels of figure |2] From 
the expressions above and figure |21 we note the following points in this case: 

1. In the small 0i3 regime (tan^ 0i3 ^ 10^^), the probabilities are again independent of 0i3 but they are not 
independent of energy anymore, as the earth matter effect introduces some energy dependence. Apart from 
012 the probabilities in this regime are now also functions of Am^i through the dependence on freg and freg- 
However as is clear from figure |21 and Eqs. (|21|l and H22|l the probabilities are same for both hierarchies and 
introduction of earth effect does not help in lifting the degeneracy of the probabilities between normal and 
inverted hierarchies. 

2. As can be seen clearly from figureia P^^and P//^ become energy dependent with the inclusion of earth matter 
effect but for a fixed energy there is no dependence of these probabilities on 0i3. 

3. In the intermediate 0i3 range (10^^ ^ tan^ 0i3 ^ 10^"^) for the Pe(i^e) channel with IH(NH) as the true mass 
hierarchy, the 0i3 dependence is governed by Ph- Whereas the energy dependence of the the neutrino (anti- 
neutrino) survival probabilities is due to both Ph and the earth regeneration factor fregifreg)- Since the 
amount of earth regeneration for these channels depends on Ph (cf. Eq. (|28|l ) the degree of energy dependence 
is different for different ranges of 0i3 - there being more earth matter induced energy variation at smaller 0i3 
as can be seen in figure [21 

4. For large values of 0i3 (tan'^ 0i3 ^ 10^'' ) P^^ and Pg|^are given by the same expressions as the no earth effect 
case and do not depend on energy or earth matter effects. But P//^ and depend on earth matter effect and 
energy. 



IV. SUPERNOVA NEUTRINO SIGNAL IN THE KAMLAND DETECTOR 



A. Energy integrated event rates in KamLAND 



Neutrinos from a SN can be observed in KamLAND through their charged and neutral current reactions with the 
detector scintillator material. The total number of charged current events triggered by the SN Pe Sux arriving at 
KamLAND can be expressed as, 



R^^ = rjNx j dE,,s j dE,<j{E,)R{E,,,,E,)acc{Eu 
X F,{E,)P,s{E,) + F^{E,){\ - Pee{E,)) 



(31) 



where N is the number of target nuclei/nucleons, rj is the detection efficiency, icciEu) are the corresponding CC 
cross-sections for the Ve and E^s is the measured visible energy of the emitted positron, when the true visible energy, 
i?^^ is related to the antineutrino energy i?^ as, i?^^ = Ey — Q.7% MeV. R{Eyis,E^+) is the Gaussian energy resolution 
function of the detector given by. 



R{Eyis,Eg 



exp 



E, 



2al 



(32) 



with ao{E)/¥. = l.b%/y/{E) as given by the KamLAND collaboration 0| . Fg is the unoscillated electron antineutrino 
flux produced in the SN and F^ stands for the flux of any of the other antineutrino flavors, i^^ or Vt. The expression 
for the neutrino flux of species a, arriving at the detector is given by 



Ea {Ejy ^ Ea , C?) 



El 



5.68Ti 



exp(f^) 



(33) 



where La is the total SN energy released in the neutrino species Va , Ta is the temperature of the Va gas at the neutrino 
sphere and d is the distance of the supernova from earth. In Ea. (|33|l we assume a pure Fermi-Dirac spectrum for the 
neutrinos. As discussed before in section^ the actual SN neutrino spectrum may deviate from a perfect black body. 
This deviation can be accounted for by introducing in the neutrino distribution function a "pinching factor" 77, which 
is a dimensionless parameter defined as /Zq/Tq, where fj,a and Ta are the chemical potential and the temperature of 
the neutrino species Va |39j |. The pinching factor has the effect of suppressing the high energy end of the neutrino 
spectrum ^. 

The number of CC events induced by the SN I'e in KamLAND is given by the same expression Eq.J^U, but with 
the antineutrino fluxes Eg and F^ replaced by the neutrino fluxes Eg and F^, Pes replaced by Pee and with the 
corresponding CC interaction cross-section for the neutrinos. The expression for the neutral current events in the 
detector is given by 



R 



NC 



= T^N / dE, 



<JNc{E.) Ee{E,) + 2F^{E,) + dNc{Eu)[Fe{E,) + 2Fs{E,) 



(34) 



where ctnc{E^) and (7nc{E^) are the NC cross-sections corresponding to neutrinos and antineutrinos respectively 
and all other variable are as defined before. 



B. Detection reactions in KamLAND 



KamLAND is an isoparrafin based liquid scintillator (80% paraffin oil and 20% pseudodocumene) detector with a 
typical composition of C„i?2n- The main CC detection reactions in KamLAND are: 



Consequences of a a small deviation from Fcrmi-Dirac in the high energy tail, for neutrino oscillation parameters, have been discussed 
for instance in |20| . 



1. The CC capture of De on free protons: 

i^e + P ^ + n (35) 

The threshold energy for this reaction is 1.8 MeV. The released positron annihilates with an ambient electron 
to produce two gamma rays. The final state neutron thermalises in about 180 jisec to produce a 2.2 MeV 7 ray. 
The delayed coincidence of these 7s gives an event grossly free from backgrounds. ^ . 



2. The CC capture of De on ^^C: 



+i2c^^2B + e+ (36) 
c + e- + i7e (37) 

The threshold antineutrino energy for the capture on ^^C is 14.39 MeV. The (3 unstable ^^B eventually decays 
with a half life of T1/2 = 20.20 ms. This signal is also recorded by the delayed coincidence of the 7s produced 
from annihilation of e+ produced in the final state of the first reaction, followed by the annihilation of e~ 
produced in the final state of the second reaction. 

3. The CC capture of i/^ on ^"^C: 

+12 C ^ + (38) 
i27v^i2 (39) 

The threshold neutrino energy for capture on ^"^C is 17.34 MeV and the /3 unstable ^"^N created in this process 
decays with a half life of T1/2 = 11.00 ms. The signal is again recorded by the delayed coincidence of the two 7s 
produced from annihilation of the e"*" and e~ respectively, produced in the final states of the two above reactions. 

The main NC reactions by which KamLAND can observe SN neutrinos are: 

1. The neutrino-proton elastic scattering: 

v^-p^v + -p (40) 

The scattered protons have kinetic energies of ^ MeV and hence cannot be observed in the high-threshold 
Cerenkov detectors. However in a scintillator detector like KamLAND, the energy deposition due to the protons 
can be observed in spite of the quenching of the proton energy due to ionisation 33J- The cross-section for 
i/g scattering is same as that of scattering by v^^ while the cross-section for j>e scattering is same as that of 
scattering by [s^. But the scattering cross-section of the particles is different from the scattering cross- 
section of the anti-particles. However under the approximation of E^/Mp « 1, where Mp is the proton mass, 
even the particle and anti-particle cross-sections can be taken as almost equal. 



2. Neutral Current break-up of ^^C: 



C -^^^ C* + u' {y') (41) 
i2fj* ^12 (J _^ ^(i5.llAfey) (42) 



The threshold energy for this reaction is 15.11 MeV. The detection of the 15.11 MeV 7 released by the de- 
excitation of the ^2(7* gives an unambiguous signal of the supernova neutrinos, virtually free from all back- 
grounds. 

All the above reactions have distinctive signatures, allowing one to count the event rates separately. The CC 
reactions can also measure the energy spectrum of the final state electrons/positrons. In the v — p NC scattering 
reaction, the scattered proton carries the energy information and hence this reaction also has sensitivity to the 
spectrum of the incoming neutrinos. However for the NC break-up of ^^C , the emitted photon has no information 
about the neutrino energy and this is the only reaction in KamLAND with no energy sensitivity. We take the S'e + P 
absorption cross-section from j^, the v—p elastic scattering cross-section from Issj and the ^"^C reaction cross-sections 
fromtSI. 



^ This is the reaction used by KamLAND to detect the reactor antineutrinos 
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FIG. 3: Plot of the event rates with la statistical error bands in the different CC reaction in KamLAND as a function of time 
of occurrence of the Supernova for different values of 613. The curves are drawn with Am2i and 612 fixed at the best-fit value 
from the global solar and KamLAND data. See text for details. 



Finally a note on the expected backgrounds in KamLAND. The Kamioka mine has a rock overburden of 2700 m.w.e 
producing 0.34 Hz of cosmic-ray muons. Since the SN neutrino signal lasts for about 10 s, the signal for the iJ^ + p 
events is ~ 300/10 — 30 Hz, while that for the Carbon events are 10/10 = 1 Hz. Thus the signal is greater than 
the cosmic ray backgrounds. Since the CC reactions are observed through delayed coincidence techniques, therefore 
they are grossly free of background problems. The only problem may come from lovif-energy radioactive backgrounds 
for the V — p scattering reaction. However Borexino has already demonstrated that it is possible to understand this 
background. If the KamLAND collaboration want to pursue the low-energy solar neutrino program then they also 
have to understand this background properly. However it should be noted that since the SN neutrino signals last for 
about 10 s and hence is a time varying signal the background for these is almost 3 times lower than the solar neutrino 
background [ssj . 



C. Supernova neutrino signal in KamLAND 



The total number of CC and NC event rates in KamLAND for the different reactions listed in the previous subsection, 
can be obtained from the Eqs. %M\ and (|34|l . In Table we present the total number of events expected in 1 kton of 
detector mass with 100% assumed efficiency, induced by a galactic supernova at a distance of 10 kpc and located at 
the galactic centre with 5 s = —28.9° (see Appendix). The lowest row in Table Ogives the number of events for the 
case of no oscillations. Also given are the changed number of events in presence of oscillation, for both normal and 
inverted mass hierarchy, and for three different representative values of ^13, chosen from the three interesting ranges 
discussed in the previous section. For events with oscillations, we present results obtained both with and without 
earth matter effects. The earth matter effects depend on the trajectory of the SN neutrino beam through the earth 
and hence on the direction of the SN from the detector. The direction of the SN with respect to the detector depends 
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TABLE I: The number of events for the different reactions in KamLAND in presence and absence of oscillation. For the case 
with oscillation we present the events for three different values of 6*13 and for both normal and inverted hierarchy. The NC 
events are unaffected by oscillation. 

on 5s^ as well as the time during the day when the supernova explodes. In Table the number of events presented 
with earth matter effect correspond to a typical time instant for which the neutrinos cross both the earth's mantle and 
core. Figure 13 shows the variation of the SN signal in the detector as a function of the time at which the supernova 
neutrinos arrive. The number of events for the CC reactions in KamLAND are shown in figureElfor the three typical 
values of 6*13. For tan^ 6*13 = 10~^ the and Vf. events correspond to both Normal and Inverted hierarchy. As we 
increase ^13 the number of Vgivf,) events for NH(IH) does not change and continue to be given by the same band. The 
events for higher ^13 values correspond to the case of inverted hierarchy for the case of v^. events and normal hierarchy 
for the case of the CC capture on ^^C. We refer the reader to Appendix A for the details of the geometry of the 
neutrino trajectory inside the earth, which is used for the calculation of the number of events in the presence of earth 
matter effects throughout this paper. 

The no oscillation number of events for the Ve+P absorption is 383 events while the v—p elastic scattering gives 676 
events. The cross-section for the v — p NC scattering reaction is about 4 times lower than the Vc.+p CC absorption 
reaction. But since all the six neutrino species take part in the former and since the u^^ and Vt are more energetic 
than Pe, one gets a larger number of — p elastic scattering events for the case of no oscillation. For the u — p events 
presented in Table HI we have ignored energy threshold effects. If we assume a threshold of 0.28 MeV, determined 
mainly by the radioactive backgrounds, then the number of no oscillation events for the v — p elastic scattering 
reduces to 270. The energy threshold for the ^^C reactions are very high and so there are very few events in absence 
of oscillations. 

Neutrino flavor oscillations convert higher energy v^jvTiy^jvT) into i^e(i^e), resulting in an enhancement of the 
charged current events. The Ve+P CC events increase by a factor of about 1.2 - 1.7, depending on the value of 
013 and the neutrino mass hierarchy. The number of ^^C CC events increase dramatically with the introduction of 
neutrino flavor conversion. The v,, — C events increase by a factor of ~ 17 — 30, while —^^ C events could see 
an increase of ^ 2 — 4 times, over the no oscillation expected rate. From the Table |I| we see that the extent of the 
increase in the number of events for a particular reaction channel depends on the value of dn and the neutrino mass 
hierarchy. However, for a given ^13 and hierarchy, the difference between the degree of enhancement as a result of 
oscillations for different reactions can be understood in terms of: (i) relative change in the average energy of the 
v^jvf, flux due to flavor oscillations and (ii) degree of dependence of the CC cross-section on the neutrino energy. The 
relative change in the average energy due to oscillations for v,, is much larger than that for j>e. This explains why the 
V(. events show a larger enhancement than the events. Between the Ve events, we note that the relative increase for 
i/g — C reaction is more than for Vf, +p. This is because the cross-section for the v^—^'^C has a sharper dependence 
on energy than the Vf,~\-p cross-section. The neutral current event rate is invariant under neutrino oscillations since 
the NC reactions are flavor blind. 



V. PROBING 6I13 AND MASS HIERARCHY 
A. Using total event rates 

From the Tableland figure 01 we observe that the increase in number of the Vf. events for normal hierarchy and the 
j>e events for inverted hierarchy depend on the value of ^13. Therefore the total number of events observed may be 
used to probe ^13 and hierarchy. 

The figure 0] plots the total number of + p and the ^^C CC events observed in KamLAND as a function of 
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FIG. 4: Plot of the CC event rates as a function of 613 for NH and IH for KamLAND. The Icr statistical error band is also 
presented. See text for details. 



tan^ 6*13 for the case with no earth matter effect (panels 1, 3 and 5), as well as for SN neutrinos passing through 
earths mantle and core (panels 2, 4 and 6). We also plot the statistical error bands (la) assuming ^/N errors. These 
plots show that the total number of charged current events can separate between different mass hierarchies (at Icr) 
for tan^^is > 2.3 x 10^^, 3.5 x 10^^, 7.4 x 10^^ for the i?e + p, C and C reactions respectively for the 

case of no earth matter effect and for tan^ ^13 > 2.9 x 10^^, 4.4 x 10^^, 3.4 x lO"'^ for the case when we include 
earth matter effect. We notice that in the case of v^. events, the difference between the number of events for NH 
and IH is enhanced with the introduction of earth matter effect if tan^ ^13 ^ a few times 10~^. For neutrinos 
events with IH, P//^ ~ sin^ 6'i2 ~ 0.3 (cf. figure [21) for the case of no earth effect whereas with earth effect it is 
P/g^ ~ P®^ ~ {freg + siu^ ^12). We havc checked that the factor freg defined in Eq. (^i^ is mostly positive for all 
energies. This implies that P/g^ is higher for neutrinos crossing the earth so that there is less flavour conversion and 
hence from Eq. |^we see that the number of events will be lower as seen in figure 01 For NH on the other hand, 
the earth regeneration effect is suppressed by the sharp fall of Ph with (?i3 and one does not observe a pronounced 
change in the number of events due to earth effect in this range of tan^0i3. The net effect is an enhancement in 
the difference between the event rates for the two hierarchies in the case of Vf. events resulting in a better sensitivity 
towards separating between these. On the other hand, for Pg events, the difference between the hierarchies is seen to 
reduce for neutrinos crossing the earth. 

To illustrate how the total number of events can be used to put bounds on ^13 and separate between the two 
hierarchies, we draw horizontal projections on the number of events axes from the points at which the slope of the 
rate-curve changes. This is shown in figure^ by the horizontal lines marked by small English alphabets a, b, c, etc. 
Next we drop vertical projections on the tan^ ^13 axis and find out the corresponding range of ^13. The 6'i3 bounds 
are marked by the vertical projections A, B, C etc. In Table ITll we give the observed number of events between the 
different horizontal lines and the inferences on hierarchy and ^13 that one can draw from these observations for the 
three CC reactions. For purpose of illustration we give this Table for the case when the neutrinos are not passing 
through the earth. However similar conclusions can also be drawn for the case when the time of occurrence of the SN 



Reactions 


Range of observed events 


conclusions drawn 




N< 485(a) 


NH with no bound on 6*13 
OR IH with t&n^OvA < 2.3x10"^ (B) 


485(a)<N<501(b) 


NH with tan^ On > 6.2 x 10-^(0) 
OR IH with tan^ ^13 < 3.4 x 10"^ (C) 


501(b)<N<530(c) 


IH with 8.7 X W-^{A) < tan^ 6*13 < 4.7 x lO-^(D) 


530(c)<N<615(d) 


IH with 2.3 X 10-^(5) < tan^ ei3 < 2.9 x 10-^{F) 
OR IH with tan^ 613 > 8.4 x 10"^ (ii") 


N>615(d) 


IH with tan^ 6I13 > 8.6 x 10"^ (E) 




N< 29(a) 


NH with no bound on 6113 
OR IH with tan^e»i3 < 3.3xlO"^(B) 


29(a)<N<33(b) 


NH with tan^ 613 > 5.6 x 10"^ (G) 
OR IH with tan^ 613 < 4.9 x 10"^ (C) 


33(b)<N<41(c) 


IH with 6.8 X 10~^{A) < tan^ 6*13 < 1.4 x 10"*(D) 


41(c)<N<44(d) 


IH with 3.3 X 10"^(B) < tan^ 6I13 < 2.3 x 10"''(F) 
OR IH with tan^ 613 > 1.2 x W^^{H) 


N>44(d) 


IH with tan^ei3 > 4.2 x 10~^(E) 




N< 63(a) 


IH with no bound on ^13 
OR NH with tan^6ii3 < 7.4xlO"^(B) 


63(a)<N<66(b) 


IH with tan^6li3 > 1.0 x 10"^(G) 
OR NH with tan^ 613 < 1.1 x lO-^(C) 


66(b)<N<70(c) 


NH with 7.8 x 10"^ (^2) < tan^ 6^13 < 2.5 x 10"''(F) 
OR NH with tan^ 6I13 > 1.3 x 10''^ {H) 


N>70(c) 


NH with tan^ 613 > 2.5 x 10"'* (E) 



TABLE II: The inferences on hierarchy and allowed range of 613 that can be drawn from figureEJfrom the observation of total 
number of events. We present the inferences for the panels 1,3 and 5 in figure |1] a, b, c, d, e correspond to the points shown 
in figure m The bounds on ^13 correspond to vertical projections from a, b, c, d, e on the ^13 axis from the points where the 
horizontal lines meet the rate-curves. 



makes the neutrinos traverse through the earth. We see from the Table that if we consider one particular reaction 
then the number of events lying between the different horizontal lines can: 

(1) identify the hierarchy uniquely and indicate the ^13 range, 

(2) allow both hierarchies and put some bound on 6*13 for either or both of these. 

In the case when a unique solution is not obtained from one reaction it may help to take recourse to a second reaction 
channel. For instance considering the Pe+P reaction we find from Table UTI that if the total event rate, TV < 485 then 
there are two possibilities: 

(i) the hierarchy is normal and no bound can be put on ^13, 

(ii) the hierarchy is inverted and tan^ 6*13 ^ 2.3 x 10""^. 

The ambiguity in hierarchy can however be resolved if in the ly^ C reaction one observes TV > 66. For this case 
one can definitely say that the hierarchy is normal and obtain a lower bound of tan^ 6*13 > 7.8 x 10^^. With the 
identification of the hierarchy as normal, the lower bound can be improved to tan^ 613 > 5.6 x 10^^ if in the De +^'^ C 
reaction one observes 29 — 33 events. Even when the number of events observed in one particular reaction gives a 
unique solution for the hierarchy and the allowed range of ^13, it may be possible to use the observed number of events 
in another reaction effectively to further constrain the allowed 613 range. 

We would like to point out here that this is an exploratory analysis which studies the sensitivity of the total number 
of events to 6113 and the mass hierarchy. The exact values of the bounds quoted in Table ITU arc not so important. An 
appropriate analysis procedure including all the uncertainties and parameter correlations will allow to study the 
conclusions obtained in more detail and with a greater statistical significance. 



energy 


Minimum detector volume (in kiloton) 


bins 


required for different reactions 


in MeV 




l^e C CC 


C CC 


21-26 


very large 


very large 


8 


26-31 


29 


very large 


11 


31-36 


51 


5 


78 


36-41 


10 


88 


very large 


41-46 


2 


9 


5 



TABLE III: Minimum volume required for a scintillator detector to statistically distinguish the number of events between two 
scenarios: neutrinos undergoing no earth matter effects and neutrinos passing through earth's core and mantle. 




FIG. 5: The CC spectra due to +p, i>e +^ C and -|-^ C events in KamLAND. The solid curves are for the case with no 
earth matter effect and the dashed curves are for the case when the neutrinos pass through earths mantle and core. The dotted 
lines correspond to the unoscillated spectrum. We plot the spectrum for three sample values of For the v^(y^ events the 
spectra for NII(III) are identical to the one presented for tan^ &13 = 10~^ and hence these are not drawn for the higher &13 
values. 
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KamLAND detector. 
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B. Using the charged current spectrum 

For the CC absorption of De on free protons and and i^e reactions on ^^C, the detectors can observe the distribution 
of the event rates with energy of the emitted electron or positron. As the average energies of the v^/vriVfi/vT) are 
greater than the average energy of Ve{ve): the oscihated Ve{ve) will have a harder spectrum. This hardening of the 
energy spectrum of the neutrinos is manifested in the observed positron/electron events in the detector, since one has 
more high energy events compared to what would be expected in absence of oscillations. In addition, if the neutrino 
survival probabilities are energy dependent, then this would introduce a further distortion in the shape of the resultant 
spectrum. There are two possible sources of energy dependence in the survival probabilities: 

(i) matter effect in the supernova - this inducts an energy dependence in Pg|^and -Pj^^ through Ph, in the intermediate 
6*13 range; 

(ii) matter effect inside the earth - this leads to Am|]^ driven energy dependent frequency modulations. 

In figure El we show the un-binned positron/electron spectrum for the (?) Ve + p, (h) C and {iii) C 

reactions in KamLAND. The dot-dashed curve gives the spectra for the case of no oscillation, the solid curves are 
for the case where we have oscillation inside the SN but with no earth matter effects, while the dashed curves 
show the spectrum expected when the supernova neutrinos cross the earth's mantle and core. For the case where 
tan^ 013 = 10~^, the CC spectra for both and 9^ channel are identical for IH and NH (panels 1,4 and 7). For Ve{vg) 
flux with NH(IH) as the true hierarchy, the event spectra are not expected to change with ^13 in the range limited 
by the CHOOZ constraints. Therefore we present the Ve{ve) spectra for NH(IH) only in the panels 1, 4 and 7. On 
the other hand the i>e{ve) spectra for III(NII) will depend on the value of 6*13 both for with and without earth matter 
effects. In the remaining panels we show the spectra for these cases for two different values of tan^ 6'i3 - 5 x 10~^ and 
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FIG. 7; Energy histogram plots of number of CC events due to Ue+p, i^e+^^C and Ve+^^C in bins of width 5 MeV for a 50 
kilo ton scintillator detector like LENA. 



lO^'^. We also note from this figure the fast variation of the spectrum shape due to earth matter effects. 

In figure El we show the expected energy spectrum for the three possible CC reactions in KamLAND in bins of 5 
MeV width. The left-hand panels show the expected spectra without earth matter effect while the right-hand panels 
display what we expect in case the neutrino cross the earth. For both the and Ve events, the solid curve denotes 
the case of tan^ 6*13 = 10^^ and it corresponds to both NH and IH. For De{i^e) events with NH(IH) the spectra remains 
the same as that shown by the solid lines for all values of ^13, however for IH(NH) the spectra will change with 613. 
These are shown by the dashed and dot-dashed lines for two different values of tan^ ^13. The dotted lines show the 
no-oscillation spectra. 



1. No earth matter effect 

The comparison of the curves corresponding to the case of no oscillation in figures and El to the ones with 
oscillations, show the effect of oscillations in general on the expected spectrum in KamLAND. 

For a given hierarchy, say IH for (panels 1 and 3 of figure jSJ or NH for z/g (panel 5 of figure jBJ, if we compare the 
spectra for different values of tan^ 6*13, we note an increase in the number of events in higher energy bins, with the 
increase in the value of tan^ ^13. This is due to the change in the relevant survival probability due to matter effect in 
the SN. For tan^ 613 ~ 10^^, Ph — 1 and we have only the lower resonant conversion in the SN. As tan^ 6'i3 increases, 
Ph decreases and eventually becomes zero at tan^ 013 ^ 10^'^. A small Ph implies a lower survival probability, 
more flavor conversion and hence harder spectrum for the incoming (anti)neutrinos. This change in the shape of 
the spectrum with ^13 can be used to limit the value of tan^ ^13 for a given hierarchy. Also, since for NH(IH), the 



Vei^e) spectrum are independent of tan^ 613 and correspond to the the curves shown by bold sohd hues in figure|B| a 
comparison of the observed spectrum against the expected ones can shed some Hght on the neutrino mass hierarchy, 
if true value of tan^ 6*13 is ^ 10~^. Note however that the value of tan^ 6'i3 and the neutrino mass hierarchy come 
correlated. For instance, if the observed spectrum corresponds to the one shown by the solid black line in panel 1 of 
figureElfor +P events, then it would be difficult to decide whether the hierarchy is inverted and tan^ 613 ^ 10^^ or 
the hierarchy is normal with no bound on ^13, using only this reaction. However, if in addition the observed spectral 
shape corresponding to electron events generated by the z^e C reaction conforms to the one shown by, say, the 
dot-dashed line in panel 5 of figure IHl then we can say that the mass hierarchy is normal and one can put a lower 
bound on tan^ ^13. 

2. Including earth matter effeet 

Matter effects inside the earth produce Am^i driven oscillations [TsL I20I I25I |28| which we see in the dashed 
curves shown in figure and right-hand panels in figure El These high frequency oscillations in the energy spectrum 
are superimposed over the dominant spectral distortion due to transitions in the supernova. Observation of these fast 
modulations in the observed spectrum would give an unambiguous signature of the presence of earth effects, which 
can be used to constrain the hierarchy and ^13. In particular, we note that there are no earth induced wiggles in the 
resultant i?e('^e) spectrum for tan^ ^13 ^ 10~^ for IH(NH). However, if the hierarchy is normal(inverted) then even if 
tan^ 013 ^ 10""^, one would expect earth matter effects in the Pe(i'e) spectrum as is evident from panels 1, 4 and 7 of 
figurelSland the solid black lines in figure|Bl which show the spectra for these cases. Thus the detection of earth matter 
effects in the antineutrino (neutrino) channel for tan^ ^13 ^ 10~^ would clearly indicate normal(inverted) hierarchy 
for the neutrino mass spectrum. If we use the earth matter signatures from just the neutrino or the antineutrino 
channel alone, we would need some prior information that tan^ ^13 ^ 10~^. This could come either from the total 
rates signature of the supernova signal itself (as discussed in the previous sub-section) or from some other terrestrial 
experiment. However, if we use the neutrino and antineutrino channel simultaneously and look for earth matter 
effects, we could pin down both the neutrino mass hierarchy and the range of tan^ ^13 ^ 10~^, if this indeed was the 
true range of tan^ Oi^. For example, if we observe no earth matter dependent modulations in the neutrino channel 
and at the same time observe the earth matter generated wiggles in the antineutrino channel, we could surely say 
that the hierarchy is normal and tan^ ^13 ^ 10"'^. 

3. Detection of earth matter effect induced modulations 

The figure 13 reflects the fact that the earth effect induced modulations are more rapid but of less amplitude in the 
low energy regime, whereas in the high energy domain they are less rapid but of greater amplitude. Since in the low 
energy regime the modulations are very frequent in energy, they are averaged out even in a very narrow energy bin. 
So it is not possible for any detector to observe the "earth matter effect induced modulations" in the lower energy 
bins (cf. figure EJ. On the other hand, towards the higher energy regime, the modulations are highly peaked and 
more wide (roughly the same order as the width of the energy bins we take), and hence may be observed in the event 
distribution, as can be seen in figure If the energy resolution is poor which means a measured electron/positron 
energy corresponds to a wide range of true energy (and hence neutrino energy) then even the wide modulations 
at high energy domain may be averaged out within this range and consequently these modulations can't be seen 
even in the higher energy bins. So a detector needs a very good energy resolution to detect the earth matter effect 
induced modulations. KamLAND being a scintillation detector fulfills this criterion very well and we note from the 
right-hand panels of figure [HI that KamLAND does see the earth induced modulations in the energy spectrum at the 
high energy end. However, even though the event rate in KamLAND for the D,, + P reaction is not very low, the 
statistics in KamLAND in the individual spectral bins for the ^^C reactions are poor with the present volume of 1 
kton. This could therefore preclude any statistically significant signature for earth matter effects in KamLAND at 
least in the ^^C reaction channels. We therefore need larger KamLAND like scintillation detectors for observing the 
earth matter effects in the detected supernova spectrum, where these ^^C reactions could be probed to look for earth 
matter effects. In Table iHll we give the minimum detector mass required for a KamLAND type scintillator detector 
to observe statistically significant earth matter effect induced wiggles in various energy bins. The Table IIIII shows 
that in general we need lower detector mass to detect the modulations in the higher energy bins. For the lower energy 
bins (which have a better statistics) the i^e C reaction is most promising and can detect the earth matter effect 
modulations with a much smaller volume detector. 
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FIG. 8: Plot of ratio of Ve(i'e)+ ^'^C CC event rates to {v+^'^C) NC event rates as a function of ^13 for NH and IH, including 
the statistical uncertainty 



4. The CC spectrum for a 50 kton detector 



Even though from figure El we see that the statistics for the + p reaction is not very poor in KamLAND, this 
reaction feels lesser impact of the earth generated modulations of the spectrum in the energy bins near the peak. On 
the other hand the ^^C reactions have greater fluctuations over the no earth effect spectra in the energy bins near 
the peak. Thus the ^^C reactions seem to have a better potential to probe earth effects in the resultant spectrum. 
But one needs a larger detector to study the energy spectra of these reactions and make a statistically significant 
statement. In figure [7| we show the corresponding expected energy binned spectrum in a 50 kton liquid scintillation 
detector. Recently a large liquid scintillation detector, LENA, has been proposed which would have a total mass of 50 
kton and a fiducial mass of about 30 kton j^. We use the largest possible detector mass for the sake of illustration. 
This figure is a "scaled up" version of figure [51 and all the features discussed in the context of figure [HI appear in a 
much magnified form here. In particular, the fast oscillations imprinted on the spectra ^^C corresponding to and 
Vf, reactions can be combined to lead to unambiguous determination of the hierarchy and limit on the value of ^13, 
as discussed before for KamLAND. For LENA it can be verified with a good statistical precision. The Ve+P process 
in LENA would have a statistics comparable to that expected in SK - with much better energy resolution and hence 
greater spectral power. 



VI. SN UNCERTAINTIES AND CC/NC RATIOS 

In the previous sections we have considered the effect of statistical errors only. However there are uncertainties 
emerging from our (lack of) knowledge of neutrino fluxes in a supernova due to: 
• uncertainty in the total value of neutrino luminosity 
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FIG. 9: Plot of ratio of v^iy^^ ^^C CC event rates to {u + p) scattering event rates as a function of 613 for NH and IH, with 
the statistical uncertainty. 



• uncertainty in the luminosity distribution between the 6 species 

• uncertainty in the average temperature of each neutrino species 

• uncertainty in the value of pinching factor in the neutrino spectra which leads to deviation from the Fermi-Dirac 
distribution 

In our calculation so far we have used — Lg — — where x can be fi or t. However these can vary as 
Le/e — (0.5 — 2)La;/2 with Lx and being equal [37l|. The average energies of i^e, Vf. and can vary in the range 
(7 - 18) MeV, (14 - 22) MeV and (15 - 35) MeV 'ST'f^The pinching factor can take values between (0 - 2) for v,^ and 
(0 — 3) for Vf. and v^- Finally the absolute luminosity can be uncertain in the range (1 — 5) x 10^^ ergs. Incorporation 
of these uncertainties can change the above conclusions. These uncertainties are correlated and a proper treatment 
requires a analysis incorporating these uncertainties in a correlated fashion. In this section we explore if it is 
possible to construct variables which are free from some of the uncertainties and introduce the CC/NC ratios. We 
study the capability of these ratios in separating between the normal and inverted hierarchies and also their sensitivity 
to 013. These ratios have the advantage that the uncertainty due to the total luminosity gets canceled out from the 
numerator and the denominator. 

As we have discussed in section 4, KamLAND is sensitive to two NC reactions (i) v C and (ii) v + p. In figure 
figure IHl we plot the CC/NC ratios as a function of ^13, for normal and inverted hierarchy with la statistical error 
bands. The NC reaction considered is the v C reaction. In figure |5| we plot the same thing but using the v — p 
scattering reaction as the NC reaction channel. 

It is instructive to look at the %-spread in the widths of the error bands relative to the central value for each case 
as well as for the CC rates in figure^ The % spread is defined as 

1 Rupper Rlower -, r>r> / An\ 

spread = — ^ x 100 (43) 

tvccntial 

where R denotes the values of the event rates or ratios at the appropriate limit denoted by the suffix. When we take 



Reactions 


CC events 


CC/u+^^ C NC events 


CC/u+p scatt. events 


Ue C 
fe +''^ C 


3.94(4.66) 
14.11(20.31) 
13.42(11.12) 


15.72(15.92) 
20.75(25.38) 
20.39(18.90) 


5.51(6.04) 
16.22(22.94) 
14.33(11.96) 



TABLE IV: The la % spreads defined in Eq. ^ at tan^ 613 = 10"^ for IH(NH). 



Reactions 


CC events 


CC/i/ C NC events 


CQjjv +p scatt. events 


i'^ +p 


11.81 





10.51 


V, C 


19.06 


13.32 


18.91 


V, (J 


4.73 


-1.63 


4.50 



TABLE V: The "separation" as defined in Eq. 14411 for different CC events as well as for the CC/NC ratios. The separations 
in this Table are calculated at a sample value of tan^ Q\z = 10~^. 



the ratio of CC event rates(x) to the NC event rate(y), the % spread is ^^(l/a; + Therefore larger the NC event 

rate smaller is the spread and so the v — y reactions are expected to have sma ller sp read than the v —^"^ C reaction. It 
is relevant to point out here that the % spread for the the CC rates goes as a/ and hence it will have a narrower 
% spread as far as statistical error is concerned. In Table HVl we show the % spreads for tan^ ^13 = 10"'^ for the three 
reactions. The Table shows that the % spread is least for the CC rates. However the effect of the absolute luminosity 
uncertainty (from which the ratios are free) may offset this advantage for the CC rates. 

When we study the effectiveness with which one can separate between normal and inverted hierarchies it is useful 
to define the quantity 

separation = 100 x f J^^"^ ' \ (44) 

\ ^central ' ^central / 

where for Veive) events the suffixes ii\ and in the first and the second term in the numerator will correspond to 
IH(NH) and NH(IH) respectively. The suffix "central" denote the the central values of event rates or ratios. The 
subscripts "upper" and "lower" are used to indicate respectively the (central + Icr) and (central - Icr) values. 

In Table we show this quantity for the CC rates as well as for the ratios for purposes of illustration for the 
specific case of tan^ 6*13 = 10^'^ and with no earth matter effect. It is clear from the table that at least at this value 
of 013, the CC/NC ratios w.r.t the v — j) NC scattering does an even better job than the CC rates, with the added 
advantage that the ratios are free from the uncertainties due to absolute luminosity. From figures El (and also figure 
21) we can convince ourselves that this conclusion is valid for other values of 6113 as well, if ^13 is not too small, for 
which in any case the SN neutrinos loose the sensitivity to separate between the hierarchies. As we have seen earlier 
with the introduction of earth matter effect for the Vf, events there is more earth regeneration for IH implying there 
is less conversion and the number of events for IH decreases and hence the difference between IH and NH increases 
(cf. fig. 2) as compared to the no earth matter effect case. Thus when neutrinos pass through earth, the z^e events 
gain a better capability of distinguishing between the hierarchies. 

The Table fvl and the figure (HI show that when the ratio is taken w.r.t the v — C reaction the potential of the 
the CC/NC variable to distinguish between normal and inverted hierarchies worsens for the +P and the v,, -|-^^ C 
reaction and the NH and IH becomes indistinguishable excepting the case of panel 6 in figure jSl for the v^, C 
reaction where earth matter effect helps. 

VII. SUMMARY AND CONCLUSIONS 

Although the neutrino fluxes from supernova suffer from inherent astrophysical uncertainties which can be large, 
the physics potential offered by them are quite rich. A core-collapse supernova provides an unique environment where 
one can study the propagation and matter effect of all three species of neutrinos and antineutrinos. The varying 
density in supernova allows the occurrence of two MSW resonances corresponding to both atmospheric and solar 
mass scale and if a detector is placed such that the neutrinos from the supernova traverse through the earth's mantle 
and/or core then earth matter effects also come into play. 



As far as the propagation in SN matter is concerned, the solar neutrino oscillation parameters lying in the LMA 
region govern the lower resonance and the propagation is fully adiabatic. Propagation through the higher resonance 
could be non-adiabatic depending on the hierarchy and ^13. For V,,{ve) with IH(NH) there are three possibilities 

determined by the mixing angle ^^13: 

(i) If i&v? Oiz ^ 10~^, the propagation at higher resonance is completely non-adiabatic, Ph = 1, and the survival 
probability is identical for IH and NH for both Vg and De- In this regime it is not possible to differentiate between 

NH and IH and only an upper limit on ^13 can be given. 

(ii) For lO^'' ^ tan^ ^13 ^ 10~^, the propagation is partially non-adiabatic for neutrinos (antineutrinos) and nor- 
mal(inverted) hierarchy and depends on Ph- Ph changes from 1 to depending on the value of ^13 and one can 

constrain ^13 within a range. 

(iii) For tan^ 6'i3 ^ 10""^ the propagation is completely adiabatic(P// — 0) for the J^e(j^e) and normal(inverted) hi- 
erarchy and flavour conversion is maximum. The probabilities do not vary with ^13 in the range permitted by the 
CHOOZ data. In this regime one can therefore set only a lower bound on ^13. 

(iv) The i'e(i'e) propagation in SN matter for inverted(normal) hierarchy have no dependence on Ph and the survival 
probabilities do not change with Ois. Therefore if Oi:^ is such that the value of Ph is expected to be different from 1 
for fe(j^e) and normal(inverted) then observation of total number of events can distinguish between the hierarchies in 
both neutrino and antineutrino channel. The power of discrimination between the two hierarchies is maximum in the 
region (iii) above where Ph differs maximally from 1. 

This dependence of the probabilities on ^13 and hierarchy gives SN neutrinos (at least in principle) remarkable 
sensitivity to probe very small values of ^13, much smaller than the projected limits from the reactor experiments, 
superbeams and even neutrino factories, and also the mass hierarchy. 

We study in detail the sensitivity of the the total observed event rates, in the De +P, Vf, -|-^^ C and v^. -|-^^ C reactions 
in the KamLAND detector to probe ^13 and to distinguish between normal and inverted hierarchy. We find that 
considering only the statistical errors into account, it may be possible to uniquely determine the hierarchy and give a 
bound on the allowed range of ^13. In some cases it may not be possible to disentangle the correlated information on 
hierarchy and the allowed range of ^13 from the observation of total number of events in one reaction channel. However 
a distinction may be possible using the observed event rates in another reaction. We find that it is particularly useful 
to combine the information from and fe reaction channels. 

Apart from the total event rates it is also possible to observe the CC spectrum due to (i) v^. +p, (ii) P,. C and 
(iii) Vf, -|-^^ C reactions in KamLAND. Since the average energies of muon and tau type neutrinos are more, the flavor 
converted electron type neutrinos and antineutrinos will have a harder spectrum even if the probabilities are energy 
independent. The matter effect in the SN may induce an additional energy dependence in the probabilities P^J^ 
and Pg|^ if ^13 lies in the intermediate range. As ^13 increases there is more flavor conversion giving rise to more 
neutrinos in the higher energy bins. This hardening of the spectrum can be used to put bounds on 6113. On the other 
hand events for IH and i>e events for NH do not display this extra ^13 dependent hardening of the spectrum. This 
can therefore be used to separate between the hierarchies. 

If in addition there is earth matter effect then it gives rise to fast modulations of the event rate with energy. For 
the Veiye) spectra, while these modulations are independent of ^13 for NH(IH), for IH(NH) they are decrease with 
the rise of ^13 and vanish for tan^ ^13 10~^. Therefore observation of these earth effect induced modulations give a 
handle to probe f?i3 and the true mass hierarchy. 

To observe this modulations one needs a very good energy resolution and KamLAND being a scintillator detector 
is ideal for the purpose. We also study the potential of the capture reactions of and on ^^C, for observing 
the earth matter induced fast fluctuations in the observed spectrum. However these reactions have low statistics in 
KamLAND which has only 1 kton of detector mass. We determine the minimum volume that is required to observe 
the modulations due to earth matter effect in a generic scintillation detector in various energy bins. We present for a 
planned 50 kton detector the energy spectrum for various values of ^13 and discuss how the energy spectrum can be 
used to differentiate between ^13 and also to discriminate between the hierarchies. 

We also construct the CC/NC ratios which are better as far as the SN imcertainties are concerned, since the absolute 
uncertainties can get canceled in these ratios. However since the ratios have a larger relative statistical error their 
usefulness is limited unless the statistics of the NC events are very large. The two NC reactions that the SN neutrinos 
can undergo in KamLAND are v — ^'^ C and v — p reaction. Since the latter has a far better statistics, it is superior to 
the V —^"^ C reaction for probing ^13 and hierarchy. We find that the CC/NC ratio w.r.t the v — -p scattering reaction 
even has a better power to separate between the hierarchies compared to the only CC rates. 

To conclude, in this paper we have made an in-depth study of the potential of the KamLAND detector in probing 
013 and sign of the atmospheric mass scale Ato§2 through observation of SN neutrinos. The results look promising 
and warrants a more rigorous statistical analysis simulating the supernova signal and including the astrophysical 
uncertainties with their proper correlation. 




FIG. 10: The geometry describing the trajectory of a neutrino from supernova inside the earth in reaching a detector. 



Appendix A 

For a given location of the supernova the nadir angle subtended at a particular detector also changes continually 
with time as the earth precesses about it axis of rotation with an angular speed of ||- hours" ^. The figurelTUldescribes 
this geometry. The direction of the supernova from earth can be specified by the angle 6s (the declination), which is 
the angle that the supernova direction makes with the plane of the equator. We denote the latitude and longitude of 
the detector by {a,6) and the angular speed of earth's rotation about its own axis by lo. We define an universal time 
which is counted from a moment (t=0) when the supernova is on the Greenwitch meridian. Then one can write the 
nadir angle 6'„ of the detector at time t as 

cos0„ — —{cos a cos S s cos{LLjt + S) + sin a sin 5s) (45) 

For cos6'„ < the SN neutrinos do not cross the earth, for < cos^n < 0.84 the neutrinos pass through the earth's 
mantle, while for 0.84 < cos 9n < 1 they pass through both the mantle and core of the earth. 

S.G. would like to thank Mark Vagins for useful correspondences. S.C. acknowledges discussions with Lothar 
Oberauer on the proposed LENA detector. 
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